pompExample(ou2)
good <- probe(
ou2,
probes=list(
y1.mean=probe.mean(var="y1"),
y2.mean=probe.mean(var="y2"),
y1.sd=probe.sd(var="y1"),
y2.sd=probe.sd(var="y2"),
extra=function(x)max(x["y1",])
),
nsim=500
)
summary(good)
plot(good)
bad <- probe(
ou2,
params=c(alpha.1=0.1,alpha.4=0.2,x1.0=0,x2.0=0,
alpha.2=-0.5,alpha.3=0.3,
sigma.1=3,sigma.2=-0.5,sigma.3=2,
tau=1),
probes=list(
y1.mean=probe.mean(var="y1"),
y2.mean=probe.mean(var="y2"),
y1.sd=probe.sd(var="y1"),
y2.sd=probe.sd(var="y2"),
extra=function(x)range(x["y1",])
),
nsim=500
)
summary(bad)
plot(bad)
Run the code above in your browser using DataLab